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ABSTRACT 

Future galaxy surveys hope to distinguish between the dark energy and mod- 
ified gravity scenarios for the accelerating expansion of the Universe using the 
distortion of clustering in redshift space. The aim is to model the form and size of 
the distortion to infer the rate at which large scale structure grows. We test this 
hypothesis and assess the performance of current theoretical models for the red- 
shift space distortion using large volume N-body simulations of the gravitational 
instability process. We simulate competing cosmological models which have iden- 
tical expansion histories - one is a quintessence dark energy model with a scalar 
field and the other is a modified gravity model with a time varying gravitational 
constant - and demonstrate that they do indeed produce different redshift space 
distortions. This is the first time this approach has been verified using a tech- 
nique that can follow the growth of structure at the required level of accuracy. 
Our comparisons show that theoretical models for the redshift space distortion 
based on linear perturbation theory give a surprisingly poor description of the 
simulation results. Furthermore, the application of such models can give rise to 
catastrophic systematic errors leading to incorrect interpretation of the observa- 
tions. We show that an improved model is able to extract the correct growth 
rate. Further enhancements to theoretical models of redshift space distortions, 
calibrated against simulations, are needed to fully exploit the forthcoming high 
precision clustering measurements. 

Subject headings: Methods: numerical — Cosmology: theory — dark energy 
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Introduction 



The accelerating expansion of the Universe can be explained either by a dark energy com- 
ponent or a modification to gravity. In both alternatives, the cosmic expansion history can 
be described using an effective equation of state, w(a), where a is the scale factor. If two 
models have the same w(a), then, as a consequence, it is not possible to distinguish them 
using a measurement of the expansion history alone. Structures are, however, expected to 
collapse under gravity at different rates in dark energy and modified gravity cosmologies. In 
general relativity, the growth of density perturbations depends o nly on the e xpansion his- 
tory through the Hubble parameter, H(a), or equivalently, w(a) ( jLinderl 120051 ). This is not 
the case in modified gravity theories. By using the measured expansion history to predict 
the growth rate of structure and comparing this estimate to a direct measurement, it has 
been argue d that it is possible to determ ine the physical origin of the accelerating cosmic 
expansion (ILue et al.l 12004 ; iLinderl 120051 ). If there is no discrepancy between the observed 
growth rate and the prediction assuming general relativity, this implies that a dark energy 
component is responsible for the accelerated expansion. 

Here we test this hypothesis using large N-body simulations which are the only way to 
accurately follow the growth of cosmic structure and hence to probe the limits of perturbation 
theory. Previous simulations of gravitational instability in hierarchical cosmologies have 
shown that linear theory gives a surprisingly poor description of fluc tuation growth and 



the redshift space distortion of cluster ing, even on large scales (e.g. lAngulo et al.l 12008 
Smith et al.ll2008b iJennings et al.ll2010bl ). We simulate the growth of structure in a modified 
gravity model and a dark energy model which, by construction, have the same expansion 
history. The growth rate is measured from the appearance of the power spectrum in redshift 
space. The goals of this paper are, firstly, to determine if these competing cosmologies can 
be distinguished from the distortion of clustering as measured in redshift space, using the 
simulation results, and secondly, to test theoretical models of the power spectrum in redshift 
space against the simulation results, to assess how well they can recover the growth rate. 

This letter is set out as follows. In Section [2] we review the growth of perturbations and 
describe the modified gravity model. Clustering in redshift space is measured in Section 3, 
and theoretical models are applied to describe the simulation results. In Section 4 we present 
our conclusions. 
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2. The cosmological models and simulations 

Here we recap how perturbation growth depends on the expansion history and the 
strength of gravity (Section I2.ip . before outlining the modified gravity model (Section 12. 2j) 
and our N-body simulations (Section l2.3p . 



2.1. The linear growth rate 

In the framework of general relativity (GR), the growth of a density fluctuation, 5 = 
(p(x,t) — p m )/Pm, where p m is the average matter density, depends only on the expansion 
history, H(a). Using the perturbed equations of motion, within GR, the growth of pertur- 
bations follows 

5 + 2H8-4<KG N p m 8 = 0, (1) 

where is the present gravitational constant found in laboratory experiments and a dot 
denotes a time derivative. The growth rate is / = dln<5/dlna, where 5(a) is the growing mode 
solution to Eq. [TJ Changing varia bles to q = S l a and allowing the gravitational constant to 



vary in time, denoted by G, gives (ILinderl 120051 ) 



d 2 g / ldln# 2 \ ldg 

+ 5 +o-Jl -T 1 ( 2 ) 




da 2 V 2 dlna J a da 

1 dlntf 2 _ 3 G(a) 
1 ]+ 2^a~~ 2~G^ il 

where f) m (o) is the matter density parameter. Eq. [2] shows that in GR, G(o)/Gn = 1 
and the growth of perturbations depends only on the expansion history, H(a). In modified 
gravity theories, however, the growth of perturbations depends on both H(a) and G(a). 



2.2. Time variation of Newton's constant 

Modifications to GR provide an alternative explanation to dark energy for the accelerat- 
ing cosmic expansion. Modified gravity theories can generally be divided into models which 
introduce a new scalar degree of freedom to Einstein's equations, e.g. scalar tensor or f(R) 
theories, and those which change dimensionality of space, e.g. braneworld gravity. In many 
such models, the time variation of fundamental constants, such as Newton's gravitational 
constant, Gat, is naturally present. 
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Self consistent scalar tensor theories are viable alternatives to GR and give rise to an 
accelerating expansion at late epochs. We refer to these as 'extended quintessence' models. 
Calculations which follow spatial variations in the scalar field have shown that, in practice, 
a broad r ange of these models can be e ffectively desc ribed with a time varying Newton's 
constant (jPettorino fc Baccigalupill2008l ; iLi et al.ll20!0l ). 



The var iation of Gn is constrained by various observations, such as the lifespan of stars 
( jTeller Il948l). the age of globular clusters (jdegl'Innocenti et al.lll996T). the m a ss of neutron 
stars (IThorsettl Il996l ) and the synthesis of light nuclei (jUmezu et al.l 120051 ; IClifton et al. 
20051 ). A time- varying Gn would also modify the temperature fluctuations in the cosmic 
microwave background, shifting the peaks to larger (smaller) scales on increasing (decreasing) 



Gn- This leads to a constraint on the variation of G, G/G 
f IChan fc Chull2007h . 



(-9.6 ~ 8.1) x 10 



-12 



yr" 



Here we consi der a simple model for G (IZahn fc Zaldarriagall2003t lUmezu et al.l 12005 
Chan fc Chull2007h . 



G 



H 2 G 



N 



(3) 



where 



/^o if a < a* 

i-SH 1 -^ 2 ) ifa*<a<« s 
1 if a > a s . 



(4) 



This parametrization describes a smoothly varying G which converges slo wly to its presen t 



Cui et al 



2010h . 



value, Gn, and is more physical than those based on step functions (e.g. 
The parameter, a*, denotes the scale-factor of photon decoupling and the parameters fi and 
a s quantify the deviation of G from the laboratory measured value, Gn, and the scale factor 
at which G and Gn are equal, respectively. The background evolution is given by 



H 2 



N 



+ a 



D 3^d\na'[l+w{a')} 



(5) 



Note we assume an equation of state w(a) = —1 in the modified gravity model to match 
ACDM. In Eq. El Q 

de is the ratio of the dark energy density to the critical density today. 
In the left panel of Fig. [TJ we plot the ratio of the Hubble rate for two different cosmological 
models with varying G, to the Hubble rate in a ACDM cosmology as a function of redshift. 
We chose to simulate the model with the maximum deviation of G from Gn which is still 
compatible with CMB measurements and solar system constraints (G — > G as a — > 1), which 
occurs for a stabilization redshift corresponding to a s — 1 (i.e. the green dot dashed line in 

Fig. CD. 
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Fig. 1. — Left: Ratio of the expansion rate, H(z), to that in ACDM for two modified gravity 
models specified by Eq. |U dot dashed green line: a s = 1, ju§ = 1.13 and solid black line : 
a s = 0.5, = 1.075. Right: The linear growth rate, /, as a function of redshift for ACDM 
(red solid), a modified gravity cosmology with a s = 1 and //jj = 1-13 (green dot dashed) 
and a quintessence model (blue dashed) with the same expansion history as the modified 
gravity model. The inset shows the ratio of / in the modified gravity model to that in the 
quintessence model as a function of redshift (green dot dashed line). 
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2.3. N-body simulations 



We use large volume N-body simulations to carry out the first direct test of the idea that 
dark energy and modified gravity cosmologies which, by construction, have exactly the same 
expansion history, can be distinguished by a measurement of the rate at which structure 
grows. The modified gravity model we simulate has the maximum deviation from Newton's 
constant that is compatible with observational constraints, as discussed above. We construct 
a quintessence model by fitting the expansion history to match the varying G model within 
0.25% over < z < 200. This model is con sistent with constraints on dynamical dark energy 
fcomatsu et alihood : [Sanchez et al.ll2009h . 



The simulations were carri ed out using a memory efficient version of the TreePM code 
Gadget-2, called L-Gadget-2 JSpringell 120051 ) . The simulation used N = 1024 3 ~ 1 x 10 9 
particles in a box of comoving length 1500/i _1 Mpc. The comoving softening length was 
e = 50/i _1 kpc and the present day linear rms fluctuation in spheres of radius 8 h~ l Mpc is 
cr 8 = 0.8. Simulations of extended quintessence cosmologies need to account for both the 
gravitationa l correction due to a varying G in the Poisson equation and a modified expansion 
history (see iPettorino fc Baccigalupil 120081 ) . In the modified gravity simulation, both the 
long and short-range TreePM algorithm force computations are modified to include a time- 
dependent gravitational constant. In both the dark energy and mod ified gravity simulation s 
the Hubble parameter computed by the code was also changed as in I Jennings et al.l (j2010al ). 



The linear theory power spectru m used to generate the initial conditions was obtained 
using CAMB (ILewis &: Bridle! 120021 ) . We adopt a ACDM linear theory power spectrum 
at z — 0, and use consistent linear growth factors in each cosmology to obtain the power 
spectrum amplitude at z — 200. In principle, as the quintessence cosmology could be 
classed as an early dark energy model, the linear theory spectrum should be modified in 
shape. However, as we have shown, such a change has a negligible impact on the nonlinear 
spectrum and on the ratio of the quarupole to monopole monents (Jennings et al. 2010a,b). 

To obtain errors on our measurements we ran 10 lower resolution simulations with 512 3 
particles, also in a box of comoving length 1500/i -1 Mpc, with different realizations of the 
density field. The power spectrum was computed using the cloud in cell (CIC) assignment 
scheme and performing a fast Fourier transform. For the initial conditions the linear growth 
rate for each model and ACDM was obtained by solving Eq. [2] numerically and is plotted 
in the right hand panel of Fig. [T] as a function of redshift. For all the models we used the 
following cosmological parameters: Q m = 0.26, Qbe = 0.74, flh = 0.044, hp = if /100km 
s _1 Mpc _1 = 0.715 and a spectral index of ris = 0.96 (jSanchez et al.ll2009l ). We have ver- 
ified that our modifications to Gadget-2 are accurate by checking that the growth of the 
fundamental mode in the simulations agrees with the linear theory predictions. 
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Fig. 2. — The ratio of the quadrupole and monopole moments of the redshift space power 
spectrum, P2/P0, as a function of wavenumber at redshifts z = 0, 0.5 and 1 (in order of 
ascending amplitude). The points show the N-body results, for the modified gravity model 
(green triangles) and the quintessence model (blue circles). The shading indicates the error 
on the ratio, estimated from the scatter over 10 lower resolution simulations. The horizontal 
lines show the predictions of linear theory model, with the colours having the same meaning 
as those used for the points. 
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3. Results 

We now briefly recap the models used to describe the redshift space distortion of the 
matter power spectrum and then (S 13. 21) fit these models to the moments of the power spec- 
trum measured in our simulations. 



3.1. Redshift space distortions 

The matter power spectrum in redshift space can be decomposed into multipole moments 
using Legendre polynomials. The ratio of the quadrupole and monopole moments of the 
matter power spectrum is plotted in Fig. [2j We model redshift space distortions in the distant 
observer approximation by perturbing the particle positions down one of the cartesian axes, 
using the suitably scaled component of the peculiar velocity. The simulation results show 
that this ratio has a strong dependen ce on wavenumb er. This can be contrasted with the 



linear perturbation theory prediction ( jCole et al.lll994h 



P 2 (k) 4/3/3 + 4/3 2 /7 

P {k) 1 + 2/3/3 + /3 2 /5 ' 1 ' 

where ft = f /b and b is the linear bias, which is unity for dark matter; Eq. [6] is independent 
of scale (horizontal lines in Fig. [2]). We note that, by considering redshift space distortions 
in the clustering of the dark matter, we are testing theoretical models against the simplest 
possible case. The distortions will inevitably be more complicated for dark matter haloes 
and galaxies, for which the bias factor b can have scale dependence (e.g. Angulo et al. 2008). 

In Fig [2] the quadrupole to monopole ratio increases in amplitude with redshift, due to 
the evolution in the matter density parameter. At z = there is a 2.5% difference between 
the linear theory growth rates in the two models. However, at this level, the measured ratios 
P2/P0 are indistinguishable on the very largest scales k < 0.02/i/Mpc (green dotted and blue 
dashed horizontal lines). At z — 0.5 and z = 1 the linear theory predictions for the growth 
rates in the two models differ by 4% and 6% respectively. The error on this ratio measured 
from the lower resolution simulations is shown by the shaded region in Fig. [2j 

In addition to the linear the ory model we consider two variants. The first is the Gaussian 



model (IPeacock & Doddslll994h 



P s {k, fi) = P r {k){l + /3^ 2 ) V- fcV ^ } , (7) 



where a p is the pairwise velocity dispersion along the line of sight, which is treated as a 
parameter to be fitted. We refer to Eq. [7] as the "linear theory plus damping" model. The 
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Fig. 3. — Left: The ratio of the monopole of the redshift space power spectrum to the real 
space P(k) at z = 0.5, as a function of wavenumber. Right: The ratio of the quadrupole 
to monopole moment of the redshift space P{k). The quasi-linear plus damping model is 
plotted using / = / true (/ = 1.05/ tnie ) as a red dot dashed (blue dashed) line. The lower 
panels show the ratio of Pq/Pt (right: P^/P®) using / = 1.05/ true to the same model using 
/ = /true (blue dashed line). 
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Fig. 4. — The ratio P2/P0 in the modified gravity cosmology measured from the high resolu- 
tion simulation (green points), together with three models for P2/P0, using the correct linear 
growth rate, / = /true (left), and the value of / obtained in a x 2 fit over the wavenumber 
range 0.01 < k(h/Mpc) < 0.25, / = / FIT (right). The shaded region shows the propagated 
errors from ten lower resolution simulations. The models plotted are indicated by the key: 
linear theory - blue dotted line, linear theory plus damping - black dashed line and quasi- 
linear plus damping - red dashed line. In the left panel the best fit value for a p (a v ) obtained 
in the range 0.01 < fc(/i/Mpc) < 0.25,with fixed /, was used for the linear theory plus 
damping (quasi-linear plus damping) model. 
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Fig. 5. — Measurements of the growth rate /. The results are plotted as a function of 
the maximum wavenumber used in the fit, fc majE (/t/Mpc). The symbols show the results of 
fitting to -P2/-P0 at z — 0.5 using different models: linear theory - squares, linear theory plus 
damping - circles, quasi-linear plus damping - triangles. The symbols are filled on scales 
where the model is a good description of the measured ratio. The error bars represent the 
la uncertainty. Left: we fit to the modified gravity model and aim to recover the true growth 
factor shown by the thick green horizontal line. Right: we fit to the quintessence model, 
with the target growth factor shown by the thick blue dashed line. 
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damping introduces a scale dependence into the ratio P^/Pq. The second variant model 
takes into account departures from linear theory, as well as including small scale damping 
(jScoccimarrol 12004 ) : 



P s (k,fx)= (P 55 (k) + 2fix 2 P se (k) 

+f 2 ^P ee (k)) x e -( fk ^ 2 



(8) 



where a v is the ID linear velocity dispersion and Pee and Pgg are the velocity divergence auto 
and cr oss power spectrum respectively measured from the simulations (see also lJennings et al. 
2010bl ). We refer to Eq. [H]as the "quasi-linear plus damping" model. We note that P^/Pq is 
more sensitive to changes in / than the ratio of the monopole moment of the redshift space 
to real space P{k), see Fig. [3l and, as a result, the 1 a errors for / are smaller when fitting 
to P 2 /P . 



3.2. Measuring the growth rate 

We now apply the above models to the simulation results. In Fig. m we plot the 
measured ratio P2/P0, for the modified gravity simulation at z = 0.5, together with the 
theoretical predictions. In the left panel, the correct value of / for this cosmology together 
with the best fit value for a p and a v in the range 0.01 < fc(/i/Mpc) < 0.25 was used in the 
linear theory plus damping and quasi-linear plus damping models respectively. In the right 
panel, the best fit value for / obtained by fitting over the same range of wavenumbers is used 
for all models plotted. The value of / obtained for the linear theory model is sensitive to 
the maximum value of k used in the fit. It is clear that both the linear theory and the linear 
theory plus damping models fail to predict the correct value for /, with the best fitting values 
differing by ~ 40% and ~ 6% respectively from the true value, see Fig. [5J All the models 
plotted in the right panel in Fig. H] use the value of / recovered when fc max = 0.25/i/Mpc. 
The quasi-linear plus damping model recovers the correct value of / over this wavenumber 
range to within ~ 0.64%. 

To test these models for the redshift space power spectrum further we vary the maximum 
wavenumber, k max , used in the fit and plot the recovered growth rate as a function of fc max 
in Fig. [51 With an accurate model we would recover the correct value for the growth rate / 
and the answer would be independent of the value of fc max adopted, with the only sensitivity 
to /c max being in the error on the growth rate. Fig. [5] shows that the quasi-linear plus 
damping model comes closest to meeting this ideal. Even this model breaks down beyond 
fc max ~ 0.3/i/Mpc, which suggests that the modelling of the small scale velocity dispersion 
can be improved. Most importantly, this model recovers the correct value for / and can 
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distinguish between the two cosmologies. The models based on linear theory perform less 
well. In fact, the answer depends strongly on the maximum wavenumber used in the fit. In 
Fig. Infilled symbols are plotted for scales over which the model is a good description of the 
measured ratio (i.e. X 2 / u ~ 1> where v is the number of degrees of freedom). 



4. Conclusions 

Forthcoming galaxy redshift surveys aim to resolve fundamental questions in cosmology, 
such as the origin of the accelerating expansion. We have measured redshift space distortions 
in two simulations with different cosmologies and demonstrated that a modified gravity 
model, described by a time varying Newton's constant, and a dark energy model, which 
have identical expansion histories, have measurably different growth rates. We have tested 
models for redshift space distortions including commonly used linear theory models. We find 
that models based on linear theory fail to recover the correct value of the growth rate. A 
quasi-linear model including non-linear velocity divergence terms is far more accurate and 
allows us to distinguish between these competing cosmologies. 

Even though we consider large scales, there are i mportant departures fr om linear theory 



which can only be modelled by N-body simulation (j Jennings et al.ll2010bl ). Without such 
guidance, the application of models based on linear theory could lead to systematic errors 
of the same order as the difference in / between competing cosmologies. In this event, such 
models would give the wrong conclusion about the physics driving the cosmic acceleration. 
We find that an improved model is able to recover the correct growth factor and hence to 
tell the models apart. This model can be applied to the measured power spectrum over a 
wider range of scales than those based on linear theory, making better use of the available 
data. Our tests show that a further improvement to this model is possible. Nevertheless our 
results show that with such improved models validated against simulations, the prospects of 
distinguishing between modified gravity and dark energy using clustering measurements are 
encouraging. 
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